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Abstract 
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are made between simulations on anisotropic and isotropic lattices involving 
finite momentum hadrons. We find evidence that having an anisotropy helps 
in extracting better signals at higher momenta. Initial results for the form 
factors f+(q 2 ) and fo{q 2 ) are presented with tree- level matching of the lattice 
heavy-light currents. 
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I. INTRODUCTION 



The use of anisotropic lattices with a finer grid in the temporal direction (a t < a s ) has 
been advocated for studies of correlation functions that produce usable signals only over 
a limited time range, T signah in physical units JTJ. Even when T signa i is reasonably large 
(e.g. > lfm), if one is working on coarse isotropic lattices with lattice spacing a > 0.2 fm 
only a few data points will lie within T signa i and fitting data becomes problematic. One 
is then naturally led to anisotropic lattices in order to retain the advantages of working 
with only a small number of lattice sites in spatial directions and nevertheless have enough 
temporal resolution within T signa i. Many of the recent examples of successful employment 
of anisotropic lattices have been in simulations on very coarse spatial lattices. These include 
investigations of the glueball spectrum and NRQCD studies of heavy hybrid states [EM 
and of quarkonium fine structure ||. Anisotropic lattices have also been useful in finite 
temperature studies ||[7] and have been employed as an alternate approach to simulations 
of heavy quarks In several of these examples one turns to anisotropic lattices for two 

reasons: to be able to extract a better signal, and to avoid discretization errors coming from 
large a^M or a t E. 

In this article we investigate the extent to which anisotropic lattices could also be useful in 
simulations of finite momentum hadrons. We have in mind, for instance, heavy meson semi- 



leptonic decays such as B — > tx (p) lu [ 1Q| . In order to be able to cover the full kinematic range 
of interest to experimentalists and map out the q 2 dependence of form factors (q^ = p^ — p ,fl 
with p M = momentum of the decaying B meson and p'^ = the pion momentum) one needs 
to calculate matrix elements between hadrons with large momenta. In the B meson rest 
frame, for instance, one would ideally like to simulate pions with momenta all the way from 
p' = \p'\ = to approximately p' = Mb/ '2. Correlators for hadrons with finite momenta are 
much noisier than in the zero momentum case, this being particularly true for light hadrons 
such as the pion. Hence one is dealing with a situation where T signa i is shrinking rapidly 
as the momentum increases. Even for reasonably fine lattices (e.g. at = a s ~ 0.1/m) once 
T S i gn ai falls significantly below lfm one might consider going to anisotropic lattices. 

To make the above statements more explicit, and come up with concrete examples that 
could guide us in the future, we have studied heavy meson semi-leptonic decays on both 
isotropic and anisotropic lattices with comparable coarse spatial lattice spacings around 
a s « 0.25 — 0.29/m. Simulations of semi-leptonic decays require good control over both 
two- and three-point correlators. We use identical operators and smearing functions on 
the isotropic and anisotropic lattices. This way one creates similar signals on the initial 
timeslices and can watch how they propagate in time on the two lattices. We then compare 
the ease with which physics is extracted on the anisotropic versus the isotropic lattice. We 
find that at higher momenta anisotropic simulations are superior in providing more reliable 
signals. 

Extracting a good signal is just one of the challenges one faces in lattice studies of semilep- 
tonic decays involving hadrons with high momenta. One must, for instance, also control 
large ap discretization errors. The second objective of the present simulations was to check 
how successful highly improved quark and glue actions are in removing lattice artifacts. 



2 



Working with Symanzik improved gauge actions [TTJ, D234 light quark jn| and 0(a 2 ) im- 



proved NRQCD heavy quark []13| actions we have studied dispersion relations, speed of 
light renormalizations and also heavy meson decay constants at finite momenta. We find 
continuum-like behaviour within 5% up to about p ~ 1.2GeV. Even for p ~ 1.5GeV, which 
on our coarse lattices corresponds to a s p ~ 2, deviations are 10% or less in many cases. 
Although additional tests are clearly still called for, we are very encouraged by these find- 
ings. In the future we plan to go to lattices finer than those used in the present initial 
study. Our experience to date indicates it may not be necessary to go very much finer to 
get phenomenologically meaningful results. With highly improved actions, simulations of 
hadrons with momenta as high as p ~ 1.5 — 2.0GeV are possible once a s < 0.2/m. It is not 
necessary to have a s p < 1. We also note that the findings of the present work should carry 
over from the actions used here to other highly improved actions. For instance, the most 
highly improved light quark action with better chiral properties than the D234 action would 
be the staggered action [0,0- One could also consider using heavy clover |T3j rather than 



NRQCD heavy quarks. It will be worthwhile investigating these alternative options in the 
future. 

In the next section we introduce the highly improved gauge and quark actions used in 
our study and discuss simulation parameters. In the current exploratory study, we have 
not tried to tune quark masses very accurately. On both the isotropic and anisotropic 
lattices we work with one light quark mass slightly heavier than the strange quark. On 
the anisotropic lattice we accummulated results for two heavy quark masses, one around 
the bottom quark and the other close to the charm quark. On the isotropic lattice only 
one heavy quark mass near the bottom quark was used. Section III concentrates on two- 
point function results. We compare effective masses on isotropic and anisotropic lattices and 
discuss dispersion relations and ratios of decay constant matrix elements for heavy mesons 
with and without spatial momenta. In section IV we present results for three-point functions 
relevant for pseudoscalar — > pseudoscalar semileptonic decays. Again comparisons are made 
between signals on isotropic and anisotropic lattices. The formfactors f+(q 2 ) and fo{q 2 ) are 
extracted. 



II. GAUGE AND QUARK ACTIONS AND SIMULATION PARAMETERS 

A. Gauge Actions 

We use the standard Symanzik improved isotropic gauge action including square and 



six-link rectangular loops [11 



G 4^\3< 12 < 12 < J' y> 



with 

P„u = ^Real (Tr{U,(x)U u (x + a^U^x + a u )Ut(x)}) , (2) 
—Real (Tr^^xjU^x + a^)U v (x + 2a ll )Uj l (x + a M + a^U^x + a u )Ul(x)}^ . (3) 



Rflu N c 
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(3 = 2N c /g 2 and ul is the tadpole improvement "mq" factor, for which we use the Landau-link 



definition in this article ||16|| . For each j3, ul must be determined iteratively via simulations. 
One could also consider using perturbative expressions for ul,. On anisotropic lattices one 
can drop rectangles that extend over two links in the time direction and one has, 

^(aniso) _ n ^~~< 1 J 5 P S s' 1 Rss' 1 Rs's | 




-/5Exo^^-^|. (4) 

The variables s and s' run only over spatial directions and one must now distinguish between 
temporal and spatial Landau-link tadpole improvement factors, u t and u s . xo is the bare 
anisotropy. It differs from the true or renormalized anisotropy, 

X = a s /a t , (5) 

once quantum corrections are taken into account. Just as one fixes lattice spacings through 
some experimental input in conventional lattice calculations, when working on anisotropic 
lattices one must, in addition, determine the ratio of the spatial and temporal lattice spac- 
ings, a s /a t , via some physics requirement. This leads to the renormalized anisotropy x- 
In the absence of lattice artifacts, it should not matter which physical quantity is used to 
fix x- Conversely the dependence of x 011 how it was determined provides a measure of 
discretization errors in the lattice system. In reference WJ\ the renormalized anisotropy was 



calculated for the action ^ mso > Q f eq.(f|) using both the torelon dispersion relation and the 
sideways potential method for several values of (3 and xo- Agreement was found between 
the two determinations within 3-4%. We will be using those results in this article. 



An alternate and equivalent procedure would be to write, 



(aniso) 




1 Rs's 

12 "Hf 



(6) 

(j) = x/Xo) an d adjust rj at fixed x to satisfy some physical criterion. 



B. Light Quark Actions 



The isotropic D234 quark action is given by |1 



„ I a 6 



a 



ra 
T 



a z 12 a z 1 j \ > ■ 



6 
1 

12 



4 



-ra 



(7) 



E * {7t(V* - ^s*vi 3) ) + CW • (V - ^C 3 V (3) ) + 



am 



(VI 



(2) 



1 



12 ^) + E(vf-^c 4 vf) 



i 



The quark fields \I/ C and the dimensionless lattice fields \I/ are related through 

vp = a 3/2 *„. 



(8) 



(9) 



Definitions for tadpole improved dimensionless covariant derivatives and field strength ten- 
sors are summarized for instance in the Appendix of reference [18|. We note here just the 
relation between unimproved F^ and the 0(a 2 ) improved field strength tensors F^ u used in 
D234 actions. 



F^(x) 



u; 



-(U^(x)F^ v (x + a^)UUx) + UUx - a^)F^(x - a^U^x - a^)) - (fi <-»• v) 



+ 



1 , 1 



6 ul u\ 



2)F ilv {x) 



(10) 



The last term is needed so that factors of are correctly removed from contributions to 
UF^W and WF^U that end up being four link objects rather than six link ones. Effects 
of this term are nonnegligible on coarse lattices. 

After tadpole improving the action we set the coefficients, Cu, Ci, Co and Cp equal to their 
tree-level value of unity in our simulations. We also work with r = 1. 



The anisotropic D234 action is 

t{aniso) t ^ ■> I 1 



"T>234 



, I (!( a, o 



ra 
~2 



a t a s - =l iz 



E^{7^ + ^7-(V-ic 3 V( 3 ) 

a; t X 



4 CL^CLjs 



:n) 



(2) 



- r — 2cr^F M 

4J: dp^dy 



(12) 



where we have used the spatial lattice spacing a s to rescale the quark fields according to 
eq.(pl). Note that the renormalised anisotropy x — a s/a>t appears in the quark action. In 
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a quenched calculation it is permissible to first fix \ in the pure glue sector and use this 
value for the ratio of spatial and temporal lattice spacings in the quark action. S^il nas 
its own "speed of light" renormalization term Cq. For fixed x this coefficient must be tuned 
to ensure correct dispersion relations in fermionic correlators. In the future, especially in 
unquenched calculations, it may be simplest to work with eq.(§) and eq.( [L2|) at fixed x an d 
simultaneously and iteratively adjust rj and Co using appropriate physics criteria. Another 
possibility is to use perturbative expressions for t] and Co throughout. In the present article 
we will use nonperturbatively determined Xi u l, u s and u t from references an d a 

one-loop perturbative estimate for C from reference [pL8|| . Just as in ^234, we set C 3 , C 4 
and Cp in S^IT^ equal to unity. 



C. Heavy Quark Actions 

It suffices to write down one expression for both the isotropic and anisotropic NRQCD 



actions, the former corresponding simply to % = 1 |13|J19| . We work with dimensionless 
2-spinor fields $ = a?J 2 Q c in terms of which 



S 



NRQCD 



Hq is the nonrelativistic kinetic energy operator, 



(13) 



A( 2 ) 
2x{a s M ) 

and 5H includes relativistic and finite-lattice-spacing corrections, 

a t SH = — ci — ; — — — a ■ B 

2 X (a s M ) 

+ ° 2 z( ' 2 (V-E-E-v)-fa 1 <x-(VxE-ExV) 

8(a s M ) v 7 8{a s M ) 

e ( A(2) ) 2 hc A(4) e ( A(2) ) 2 fl5) 
4 8 X (a s M ) 3 5 24 X (a s M ) 6 16n X 2 (a s M ) 2 ' 1 ; 

All derivatives are tadpole improved and, 

A (2) = £vf, A< 4 ) = £vf (16) 

i=i j=i 

V fc = V fc - ^Vf (17) 

The dimensionless Euclidean electric and magnetic fields are, 

Ek = Fki, B k = —^eijkFij (18) 

Vfc, Vfc , j = 2,3,4 and F^ v are the same as in the light quark actions of the previous 
subsection. In our simulations we again set all q = 1. 
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D. Simulation Parameters 



Simulations have been carried out on 8 3 x 20 isotropic and 8 3 x 48 (x = 5.3) anisotropic 
lattices using ensembles of 200 configurations each. Details are summarized in Table I. The 



/3, Mo and a -1 values are taken from reference [12] for the isotropic and from reference [17 
for the anisotropic lattices, respectively. The latter reference also provided xo and x- None 
of the lattice spacing determinations are very precise. For the anisotropic lattice we have 
taken the 4.503GeV quoted in [O for aj 1 from T S-P splittings and reduced it by 22%, 
approximately the amount by which light hadron a -1 values differ from those fixed by T 



splittings in quenched calculations. This is also consistent with [17] 's finding for a -1 from 
the string tension. In heavy-light physics one expects light hadron a -1 values to be more 
appropriate. The "speed of light" renormalization coefficient Cq was estimated from the 



one-loop perturbative result for x — 5.3, C = 1 — 0.45a s fI8||, using a s ~ 0.4. 



The bare light quark mass has been adjusted so that one has the same pseudoscalar-to- 
vector ratio, P/V, on the isotropic and anisotropic lattices. We found that this required 
approximately a factor of 3/2 more BiCGstab iterations to create light quark propagators 
on the anisotropic lattice compared to on the isotropic one (actual numbers are given in Table 
I). We believe the reason is critical slowing down due to an increase in the condition number 
with a^ 1 . Although our pion is still quite heavy (~ 840MeV), on the anisotropic lattice 
we did, on a few occasions, encounter problems with "exceptional" configurations. In the 
process of accumulating an ensemble of 200 configurations, 5 had to be skipped. We will see 
below that the anisotropic pion propagators are noisier than the isotropic ones for the same 
pion mass. Other mesons such as the p and B, which use the same light quark propagators 
are not affected, however. Although this aspect of simulations on anisotropic lattices is 
worrisome and needs further study, we do not, at the moment, believe that this implies 
anisotropic lattices will never be useful for light quarks. The problem has been exaggerated 
in the present calculations since we are on very coarse spatial lattices and x = 5.3 is a 
large anisotropy. In future, more realistic, simulations we plan to use finer lattices and work 
with moderate anisotropies such as x = 2 — 3. We expect to be able to go to lighter pions 
there without encountering problems. We mention that the tadpole improvement adopted 
in eq.(|T0|) leads to a larger effective "csw" than has been used in the past, and it is known 



that large csw, such as the nonperturbative csw, leads to problems on coarse lattices. This 
phenomenon appears to set in at heavier quark masses when x > 1. There were no hints of 
problems with exceptional configurations on our coarse isotropic lattice for the light quark 
masses considered. 



III. RESULTS FROM TWO POINT FUNCTIONS 

We calculated two-point correlators for the light-light pseudoscalar and vector, the B s 
and B* and, on the anisotropic lattice, also for the D s and D* mesons. We will generically 
call these the "pion", "rho", "B meson" and "D meson" channels respectively. Current 
sinks corresponding to the temporal and spatial components of the heavy-light axial vector 
currents were also considered. We employed gauge invariant smearings, separately for the 
light and heavy quarks, of the form 
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(19) 



For light quarks we used I = 10 and c sm = 1/12. We have two smearings for the heavy 
quarks both with c sm = 1/24 but with different I values, h — 2 and l 2 = 10. We found 
that l 2 works better for zero and low momentum B mesons and l\ gives better signals for 
higher momentum correlators. We accumulated data at 7 spatial momenta, (0,0,0), (0,0,1), 
(0,1,1), (1,1,1), (0,0,2), (0,2,2) and (0,0,3) in units of 27i/aL, averaging over all equivalent 
momenta. In the following subsections we will compare effective masses, present fit results 
for a small number of states and dispersion relations, and look at ratios of matrix elements 
with current sinks. 

A. Effective Masses 

We start by comparing raw data for B meson correlators on the anisotropic and isotropic 
lattices. Figs. 1. and 2. show fractional errors of the correlators versus physical time for 
several momenta. One sees that, where there is overlap, the fluctuations of B correlators 
are identical on the two lattices. Effective mass plots are presented in Figs. 3 to 5 for the B 
meson for three representative momenta. Each figure shows results in lattice units on the left 
and results in physical units, using central values of the scales from Table I, on the right. We 
label time slices such that sources are at t — 0. By looking at the anisotropic data one finds 
plateaus starting at about t = 0.3fm - 0.4fm for all momenta. This must also be the case for 
the isotropic data since the original signals are essentially the same. However, based only 
on information from the isotropic data, this is not always obvious, especially for the higher 
momenta. There is also an unfortunate upward fluctuation at around t = 0.7fm which makes 
fitting the isotropic data more difficult. One agonizes over whether to fit starting before or 
after the hump. This illustrates some of the immediate advantages of anisotropic lattices. 
Once a plateau has set in it is more easily recognized, and one is also less sensitive to one 
or two points fluctuating up or down, since there are enough other points around. 

Effective mass plots for pions and rhos are given in Figs. 6-9 and corresponding corre- 
lator fractional errors in Figs. 10 and 11. In contrast to the B meson case, one sees that 
fluctuations are considerably larger for the pion correlators on anisotropic lattices. This 
indicates that some configurations are close to being "exceptional". The rho correlators 
show no enhanced errors relative to the isotropic lattice. Despite the larger fluctuations, it 
is not difficult to fit the anisotropic pion data and fluctuations cancel to some degree upon 
calculating physical quantities such as dispersion relations or semi-leptonic formfactors. In 
Fig. 12 we show fitted energies for the anisotropic lattice pion and rho at momentum (1,1,1) 
versus t m i n /a t from single cosh fits. The pion energies have larger errors but look otherwise 
normal. In future simulations one will have to monitor pion correlators carefully. 

Meson correlators were calculated up to the center of the lattice in the time direction. 
We find that for zero and low momenta T S i gna \ covers the entire T/2 region (1.24fm for 
the anisotropic and 2.25fm for the isotropic lattice). However as the momentum increases 
T S ig n ai starts to shrink. For the B meson this occurs only at the highest two momenta 
(0,2,2) and (0,0,3). For the pion and the rho, T signa i shrinks below ~ Ifm once one goes 
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beyond momentum (1,1,1). It then becomes much harder to get reliable information from 
the isotropic lattice. For instance, it should be obvious from Fig. 9 that trying to fit the 
isotropic data is considerably more frustrating than carrying out fits to the anisotropic data. 
Fig. 13 shows fitted energies versus t m i n /a t for the isotropic rho and pion with momentum 
(1,1,1). This should be compared with the corresponding anisotropic results in Fig. 12. 
With the isotropic data, picking t m i n is extremely tricky. In past work on isotropic lattices 
we often adopted a criterion whereby if t^ in is the smallest t min which gives an acceptable 
Q- value (> 0.1) then our preferred choice for t min in a single exponential (cosh) fit to a single 
correlator would be t^ in + 2a t [p0| , plj] . This is a fairly conservative criterion leading to larger 
statistical errors than if one chooses a smaller t min . We will call fits using this criterion 
"B-fits". In the present calculation, B-fits would dictate t min /a t = 3 or 4 for most of our 
fits to isotropic data. We will abide by this t m i n for fitting B meson correlators. However, 
for the rho and the pion we will be less conservative and use t m i n /a t = 2, which is equal 
to 0.5fm in physical units. Q- values are always acceptable already for t min /a t = 2. We call 
this latter choice "A-fits". Based only on the isotropic data, one would be hard pressed to 
argue why A-fits should be preferred over B-fits. In the present comparison study, however, 
we have additional information from the anisotropic lattice that tells us that plateaus have 
set in by t ~ 0.3 — 0.4fm, i.e. before t = 0.5fm, and hence A-fits should be fine. We should 
also note that correlated fits were used throughout our analysis. 

The results presented in this subsection demonstrate that for the pion and rho correlators 
there is a clear advantage to anisotropic simulations starting with momenta around (1,1,1). 
For B correlators benefits start around (0,2,2). Below we will present physics results ex- 
tracted from the isotropic and anisotropic lattices, concentrating on momentum dependent 
quantities. For the B meson on both lattices and the pion and rho on the anisotropic lattice 
we can go to the highest momentum (0,0,3). On the isotropic lattice, pion and rho results 
can only be extracted up to momentum (0,0,2) using A-fits. With B-fits, only momenta up 
to (1,1,1) can be reached and errors are larger than with A-fits. 

Another lesson to be drawn from the present exercise is the importance of good smearings. 
We can get away here with single exponential (cosh) fits, because our smearings are rea- 
sonable and plateaus set in well within T signa i. In the case of poorer smearings it would 
still be advantageous to be on an anisotropic lattice with many data points within T signah 
since that would facilitate multi-exponential fits. Going to lattices with longer time extent 
in physical units in search of a plateau will not work for high momentum correlators since 
Tsignai stops well before the end of the lattice. It is much more profitable to explore better 
smearings and/or increase the number of points within T signal- 



B. Dispersion Relation 



A convenient way to investigate how well lattice simulations are reproducing relativistic 



dispersion relations is to consider the following quantity [12 



C(p) 



lE 2 (p) - E 2 (0) 

p2 



(20) 
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where p = \p\ and E(p) is the total energy of the particle. In a relativistic theory C(p) = 1 
for all p. C(p) is shown for the pion and the rho in Figs. 14 and 15. One sees that on 
both the isotropic and anisotropic lattices, continuum behavior is observed to better than 
~ 5% for momenta of up to 1.2GeV at the minimum and, in the case of the rho, even 
beyond that. Hence, it appears one can simulate particles with momentum as high as 
p ~ 1.5/ a s to 2.0/ a s without introducing large discretization errors. One notices slightly 
better behavior on isotropic compared to on anisotropic lattices for the smaller momentum 
points. The anisotropic results are more sensitive to the tuning of Co, for which we have used 
the one-loop estimate. We have studied the effect of this one-loop correction by calculating 
anisotropic light quark propagators with Co = 1 on a subset of 50 configurations (at the same 
time adjusting a t mo to get the same pion mass). The fancy squares in Fig. 16 show results 
for C(p) for Co = 1, i.e. without the one-loop correction. They are compared with results 
from 50 configurations using Cq = 0.82. For completeness we also include corresponding 
points from Figs. 14 and 15, which are based on the full set of 200 configurations with 
Co = 0.82. Although effects from the one-loop correction are not large, nevertheless, they 
are crucial for getting the ~ 2% to ~ 5% agreement with relativistic behavior. Even if we 
had attempted to tune Co nonperturbatively, it would have been difficult to do much better 
simultaneously for both pions and rhos, than our perturbative estimate. The errors on C(p) 
for anisotropic pions (from a bootstrap analysis) are smaller than one might have expected 
based on the large fluctuations noted above in pion correlators. We are evidently observing 
large cancellations between fluctuations in the zero and nonzero momentum correlators. 

The isotropic results in Figs. 14 and 15 are based on A-fits, defined above. In Fig. 17 we 
compare with results using B-fits. Results from the two different fits are consistent. But, 
as expected, errors are much larger for B-fits and one cannot go beyond the three lowest 
momenta using them. Unless specified otherwise, for the rest of this article we will use 
A-fits for pions and rhos on isotropic lattices. The isotropic results in Figs. 14 and 15 are 



also consistent with findings in reference [12], where comparisons were made with clover 



quarks at similar lattice spacings. In those studies C(p), at lattice spacing a = 0.25 fm, was 
consistent with unity within errors up to about p = 1.2GeV for the D234 quark action and 
exhibited 10% ~ 20% deviations from unity for the clover quark action. 

C(p) is not applicable for heavy-light mesons involving NRQCD heavy quarks. The NRQCD 
action omits the heavy quark rest mass and meson correlators fall off with an energy E S i m (p) 
that differs from the total energy E(p). However one has, 

5E(p) = E sim {p) - E sim (0) = E(p) - E(0), (21) 

and one can define the "kinetic" mass, M^ n = M2 through, 

M km = (p 2 - 5E 2 (p))/(2 5E(p)). (22) 

If E sim (p) has the correct momentum dependence, then M kin should be independent of the 
momentum used in the RHS of eq . (|22|) . Fig. 18 shows M^n for several mesons versus the 
momentum of the correlator from which it was extracted. On the anisotropic lattice we have 
data for two heavy- light mesons, both the B s and the D s . The latter meson uses NRQCD 
charm quarks. One sees that for all mesons, results for M kin are independent of momentum 
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within errors up to about p ~ 1.5GeV. For the pion and the rho we can also compare Mki n 
with the rest mass Mi = E(0). The deviation of M/, in from Mi reflects the deviation of 
C(p) from unity in Figs. 14 and 15 and, like the latter, is very small. 



C. Decay Constant Ratios 

Heavy meson semi-leptonic decay formfactors, which will be the focus of the next section, 
require matrix elements of heavy-light currents between hadronic states with and without 
momentum. It is worthwhile considering first simpler matrix elements of currents between 
mesons and the vacuum and studying their momentum dependence. Such matrix elements 
are relevant for meson leptonic decays. Starting from the usual definition of the decay 
constant (in Euclidean space), 

(0\Ap\B,p)=p ll f B , (23) 
{pn = (iE,p)) one can form the ratio, 



(Q\Ao\PS,d/y/E(p) ^E{p) 



(0\A \PS , p= 0)/VMp~s~ ^Mp~ s ' 
This ratio was studied recently on finer lattices using a less improved action |22| (similar 



calculations were done several years ago with relativistic heavy fermions in [23]). In the 
present article we will use the following current for the temporal component of the axial 
vector current, 

A ^J^=y l5lQ Q, (25) 

where \I/ is the light quark field and the heavy quark 4-spinor Q has the NRQCD 2-spinor $ 
as the upper two components and zero for the lower two components. The superscript (0) 
signifies that J$ is the zeroth order term in an 1/M expansion for the axial vector current. 
The LHS of eq.(p4|) is replaced by, 



(0\J^\p = 0)/VMp- s 



In reference |22[ 1/M current and one-loop matching corrections were included in the ratio 
and seen to have only a small effect relative to using just the zeroth order current for 
heavy quark masses around the 6-quark mass. Matching calculations for the actions of the 
present article have not been carried out yet so we are forced to use the simple ratio (p) 
here. Figs. 19 and 20 show for the B s and D s leptonic decays, compared with the 
expected continuum behavior of the RHS of eq . (^4|) . One sees good agreement for most of 
the momentum range studied. Only at the highest momentum (> 1.5GeV) does one see 
~ 15% deviations for B s mesons. One should be able to reduce these errors below ~ 10% 
by going to slightly finer lattices. Since we do not include higher order (in p/M) current 
corrections one expects agreement with full continuum QCD behavior to be worse for the 



D s meson, especially at higher momenta [2^]. The better agreement found in Fig. 20 as 
compared to in Fig. 19 is, hence, fortuitous. 
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D. Zero Momentum Spectrum 

In the process of studying momentum dependence of meson correlators we have also 
accumulated some zero momentum spectrum results. They are summarized in Table II. 
The first errors are statistical and the second represent errors coming from uncertainties in 
a -1 which we take to be roughly 10% not including quenching effects. In this exploratory 
study we will not try to estimate other systematic errors. One sees that the "pion" and 
the "rho" masses are very close to each other on the isotropic and anisotropic lattices. The 
isotropic B s meson is about 250 MeV too heavy and the anisotropic one about 600 MeV 
too light compared to experiment. On the other hand, the uncertainty coming from a~ l is 
at the 500-600 MeV level. 

It is amusing that the heavy-light hyperfine splittings, which have been the bane of lattice 
heavy- light spectroscopy in the past, agree so well here with experiment for both the B s 
and D s mesons. One should not make too much out of this, however, until results at 
other lattice spacings have been obtained. At the moment we have no estimate of the size 
of scaling or unquenching corrections to this quantity. If one corrects approximately for 
the incorrect heavy quark mass in our simulations and multiplies by a factor of [calculated 
meson mass] /[experimental meson mass], then the entries in Table II for the B*—B s hyperfine 
splitting, are modified to 48MeV or 40MeV for isotropic and anisotropic lattices respectively. 
Similarly the D* — D s splitting becomes 148MeV. All these numbers are much larger than 
and in better agreement with experiment than in previous lattice estimates [PT] , P5| -|2"9|| . One 
difference between the present and previous calculations lies in the last term in eq.fllOD which 
enhances F^ v . It is hard to imagine, however, that this term alone can be the whole story. It 
will be interesting to see what happens after systematic errors have been investigated more 
thoroughly and various corrections to the present calculation have been incorporated. 



IV. RESULTS FROM THREE POINT FUNCTIONS 

For pseudoscalar to pseudoscalar semi-leptonic decays one is interested in matrix elements 
of the heavy-light vector current V 1 . The matrix elements are then parametrized in terms 
of form factors /+ and /_ (or /i and /o) . 

(Tr(p')\V»\B(p)) = f + {q 2 W+p'n + /-(g 2 )(p* -p**), (27) 

where q 2 = (p 11 — p' M ) 2 . We will write formulae for B decays but they apply also to D decays. 
fi and fo are defined as : 

q 2 

fo = f++ (M 2 -Mi) f - ' /l = /+ (28) 

Our main goal in this section will be to compare the quality of signals for matrix elements 
such as those in eq. (|27D between isotropic and anisotropic lattices. Since the light quarks in 
this study are still heavier than the strange quark, the physical situation we are simulating 
will be closest to B s or D s decays into Kaons. 
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A. Current Matrix Elements 



In order to extract the matrix element of eq.(B7p one starts from the following three- 
point function, G^(p,p',t B ,t), on the lattice (for technical reasons it is more convenient to 
consider the time reversed matrix element) 

G^(p,p\ t B , t) = E E e-^ s e^-^{^ B {t B , x) Vfo y) *+(0)|0>. (29) 

x y 

$J. and & B are interpolating operators used to create the pion or B meson respectively. 
is the dimensionless Euclidean space lattice heavy-light vector current. It will be defined 
more precisely below. It is related to the continuum Minkowski space through, 



C = a*V4 0) £(/^- (30) 

Z^ is the tree-level wave function renormalization for lattice light quark actions. It is dis- 
cussed for the isotropic and anisotropic D234 actions in the Appendix. £ (//) is the conversion 
factor between Euclidean and Minkowski space quark bilinear currents which is necessary 
due to the different 7-matrix conventions in the two spaces. £(0) = 1 and = —i, 
k = 1,2,3. t B denotes the time slice at which the B meson operator is inserted. In the 
simulations t B is kept fixed and we vary t, the timeslice of the current insertion, between 
and t B . Physics is extracted from those timeslices where the corresponding two-point cor- 
relators are dominated by the ground state and where e~ Ew ^ T ~ t ' can be ignored relative to 
e~ Ent , T being the time extent of the lattice. If these conditions are satisfied the three-point 
correlator (E9|) becomes, 



[2E B a i s ) {2E n ai) 

The exponential factors in eq.(^TJ) can be removed by dividing with the appropriate two- 
point functions. 

G { B \p,t) = Ec"**<0|* B (t,5)$t(0)|0> 

e - E -f»t = c, BB e~^\ (32) 



|(S(p)|^|0)p ™™ t 



(2E B a 3 s 

G®(P\ t) - [e-*-* + e~ E ^} » U e~^. (33) 

[2E 7r a a s ) 

The matrix element of the continuum current V 1 can now be obtained from, 

G^(p,p',t B ,t) 
G ( B '(p,t B -t) G^(p',t) 



T^(p,p',t B ,t) ee C(p) " ' (a) r . .J CbbC** (34) 



yf&fWpm = (35) 
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Eqs. (0) and ([35|) relate the three- and two-point functions evaluated in our simulations 
to the continuum matrix elements of eq.(^). We need to specify now the lattice current 
that enters into the three-point functions. It is defined in terms of Euclidean space 
7- matrices 7 M = £ (/i) 7(Mm& ) ■ Since we use the NRQCD formulation of heavy quarks, 
becomes an expansion in 1/M, the inverse of the heavy quark mass. After matching to 
continuum QCD one has, 

^ = £cf>jg. (36) 
j 

For j = one has the zeroth order, 0((1/M) ), current, 

J$}=^ 1 ,Q, (37) 

with the fields $ and Q defined as in eq. (p5|) . Higher order currents (j > 0) are listed in 
reference |3(J. The matching coefficients Cj V ^ are not yet known beyond tree-level so we will 
work with Cq ^ = 1 and all other C/s equal to zero. One of the 1/M current corrections 
j\P has c[ V ^ = 1 + 0(a s ) and also contributes at tree-level. However, its matrix elements 
include power law terms that will not be cancelled unless a proper one-loop calculation has 
been carried out. Hence, we do not include Jy contributions in the present study. 

We have evaluated rt 3 ' (p, p', £g, t) of eq. (|34]) for several pion momenta p' ranging from (0,0,0) 

to (0,0,3) in units of 2n/(aL). The B meson momentum was always set equal to zero. On 

the anisotropic lattice we used t_e/a t = 24 and on the isotropic lattice tsjat = 10. Figs. 
(3) 

21 and 22 show T versus t in physical units for zero momentum and for three nonzero 
momenta. One sees that, where there is overlap, anisotropic and isotropic lattices give 
consistent results. One is interested, of course, in the region where Iq is independent of 



t implying that the simple t-dependence in eq.(|3~l~D and the last expression in eq.([33|) is 
justified. can then be identified with the asymptotic matrix element (V°) of eq. (j35|) . 
Again the crucial question for the higher momentum isotropic results becomes whether one 

would believe in the presence of a plateau if one did not have the comparison anisotropic 

(3) 

data. To illustrate this point we show T versus t in lattice units separately for isotropic 
and anisotropic lattices in Figs. 23 to 26. One might still feel comfortable extracting a 
signal from the (1,1,1) isotropic data. One would be hard pressed, however, to claim that 
a plateau has been established at momentum (0,0,2) based solely on the left hand plot in 
Fig. 26. Note that T signa i, the time range over which statistical errors are under control, is 
about Ifm for momentum (1,1,1) and has shrunk to about 0.5/m for (0,0,2). We saw in the 
previous section that individual two-point correlators had reached a plateau by 0.3 ~ 0.4/m. 
If sufficient data points could be introduced between 0.3/m and 0.5/m, then one should be 
able to extract meaningful results for (V°). Hence, once again one sees an advantage to using 
anisotropic lattices starting with momentum (1,1,1). This is the same situation as with the 
pion and rho correlators described in the previous section, not a surprising finding since 
statistical errors in r^ 3 ' are dominated by the pion correlator and not by the B correlator. 

From the region where is independent of t, one can extract (V^). We show (V°) and 
(V k ) in Fig. 27 as a function of the pion momentum p' . For points at the largest momenta on 
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the anisotropic lattice one could be seeing some discretization effects. Only a more careful 
analysis involving simulations at several lattice spacings and/or further studies with nonzero 
B meson momenta will be able to shed more light on this. Here we are concentrating mainly 
on whether signals can be extracted, postponing scaling studies for the future. This is in 
contrast to the two-point correlator studies of the previous section where several continuum 
expectations based just on simple Lorentz symmetry considerations could be tested. 

B. Form Factors f+(q 2 ) and fo{q 2 ) 

From (|35|) , (^7j) and (|28"D one can extract the form factors f+(q 2 ) and foil 2 )- Isotropic 
lattice results for the B s meson are shown in Fig. 28 and anisotropic lattice results for the B a 
and D s mesons in Figs. 29 and 30. The kinematics, including the range in q 2 that is covered, 
depends on the meson masses of Table II and differs for the two types of lattices. The errors 
in the three figures are statistical and come from a simultaneous bootstrap analysis of the 
V° and V k three-point functions and the B (or D) and pion correlators. 

Since we have at present only tree-level matching of the heavy-light currents and also have 
not tuned (or extrapolated in) the heavy and light quark masses, the above form factor 
results cannot be applied yet to phenomenology. What is important, however, is that with 
just 200 configurations it was possible to obtain form factors for a nontrivial range in q 2 and 
that we were able to demonstrate the advantages of anisotropic lattices in calculations of 
this kind. By increasing statistics in the future one should be able to get good data over an 
even wider range in q 2 . At that point one would also want to go to finer lattices and work 
with improved currents so that finite momentum errors are minimized even at the largest 
momenta for which signals can be obtained. 

V. SUMMARY 

This article investigates the extent to which anisotropic lattices can help in extracting 
better signals from two- and three-point correlators involving high momentum hadrons and 
whether they can play an important role in studies of semi-leptonic heavy meson decays 
into light hadrons. To address this question we have carried out simulations of heavy meson 
semi-leptonic decays, in parallel, on isotropic and anisotropic lattices. In order to have a 
meaningful comparison, we work with similar coarse spatial lattice spacings and use identical 
sources and smearings on the two lattices. 

We find that it is considerably easier to extract reliable signals from anisotropic simulations 
once the light hadron momentum reaches (1,1,1) 27r/(a s L) or higher. This advantage may 
not be so obvious just by comparing Figs. 28 and 29 and one needs to go back to figures 
such as Fig. 22 to fully appreciate how much the anisotropy is helping here. The last 
point (at the smallest value of q 2 ) in Fig. 28 comes from the isotropic data in the left 
hand plot of Fig. 22 (see also Fig. 25). It is because the first three isotropic points in 
Fig. 22 agree with the anisotropic data, that one feels confident about the form factor 
results in Fig. 28. Without the anisotropic data, one would have to allow for a considerable 
additional systematic error, which one might call "fitting error" or "t m j„ dependence error" , 
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when presenting isotropic results. Hence, the main conclusion from the present work is 
that anisotropic lattices definitely improve signal quality and should be considered in semi- 
leptonic decay studies, especially if a large range in q 2 is of interest. 

The advantages of anisotropic lattices come at a certain price. For instance, such lattices 
require more sites in the time direction. Also light quark inversions take more iterations 
in order to get the same pion-to-rho mass ratio P/V. In the present simulations, the cost 
increase from just these two effects meant a factor of 2.4 x 1.5 = 3.6 in CPU time. Working 
with anisotropic lattices also requires tuning of more parameters. At a minimum, two addi- 
tional parameters, r] in the glue action and Co in the D234 quark action must be determined 
nonperturbatively or perturbatively. It is important that efficient procedures for carrying 
out such tunings be developed. Another drawback in the current anisotropic simulations was 
the increased susceptability to exceptional configurations. If it were not for this problem and 
the larger fluctuations in pion correlators, the anisotropic results for the form factors fo{q 2 ) 
and f+{q 2 ) would have been of even higher quality relative to those coming from isotropic 
lattices. There are several ways one could try and ameliorate this problem in the future. 
Just going to finer lattices and working with more moderate anisotropies, \ < 5.3, should 
help. One could also explore other light quark actions that have better chiral properties 
than the clover or D234 actions, such as highly improved staggered fermions domain 
wall fermions |3T| or the twisted QCD approach ||32|| . 



Ultimately one will have to weigh the additional costs associated with anisotropic lattices 
against the likelihood that an isotropic simulation with much higher statistics will allow 
us to approach the high momenta we are seeking. Based on our (and also other peoples') 
experience to date, we do not believe isotropic lattices can be competitive at high momenta. 
Signal-to-noise of high momentum correlators will be decreasing exponentially with t. When 
a t is large it will be very costly to move even one additional point from the noise into T signa i. 
In order to do so one will need to reduce errors by roug hly a factor of e at{E{p) - Ea \ where 



EO is the ground state energy [Q. In other words, one will need an increase in statistics 
by a factor of N sta t = \e at ^ E<ypS} ~ E0 ' > \ 2 . Typical numbers for the current isotropic simulation 
for momenta starting with (1,1,1) and higher, would be N stat ~ 4 — 10. This means a factor 
of 16 - 100 if one wants to go from just one or two points in T signa i to a marginally useful 
Tsignai of 3 to 4 points. This is much more than the cost of going to anisotropic lattices. 

We mention another unrelated advantage of anisotropic actions. One difference between 
<Sg so) and S { a niso) is the omission in the latter of rectangles that span two links in the time 
direction. Similarly includes higher time derivatives that are absent in S^IT^ ■ As 

a consequence, for the same amount of improvement in spatial directions, the anisotropic 
actions suffer less from lack of reflection positivity and/or presence of ghosts. In addition to 
being theoretically cleaner this means perturbation theory is more straightforward for the 
anisotropic actions. One such example is the Z^> calculation in the Appendix. 

In the course of this study we accummulated a wealth of two-point correlator data at finite 
momentum. This enabled us to compare lattice results for dispersion relations and leptonic 
decays of finite momentum heavy mesons with continuum expectations. We found that 
momentum dependent discretization errors were under control and less than 10% up to 
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about a s p m 2. Discretization errors at high momenta in three-point functions have not 
been critically assessed to date. However, assuming a situation similar to the one found 
with the two-point functions, prospects for simulating hadrons in semi-leptonic decays with 
momenta as high as 1.5 — 2GeV look promising. Only slightly finer lattices (a^ 1 > lGeV) 
than those of the present work may be required. This would have to be coupled with 
an anisotropy of \ ~ 2.5 or higher in order to be able to extract a signal at those high 
momenta. The experience gained in the present work will be indispensable when picking 
optimal simulation parameters in the future and going onto more realistic calculations. 
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APPENDIX A: TREE-LEVEL WAVE FUNCTION RENORMALIZATION FOR 

D234 QUARK ACTIONS 

In this appendix we sketch the derivation of the tree-level wave function renormalisations, 
, for the light quark actions used in this article. Results for S^u^ have already appeared 
in reference [[HJ . One starts from the requirement that the propagator for a zero momentum 



quark have the form, 

f-vr/at dp Q 



1 ^e™ t G(p o ,p = 0) 



-ir/at 2lT 

= Z q e- M ^-±^ + ... . (Al) 



G(p) is the momentum space propagator and M\ denotes the pole mass. The dots refer 
to lattice artifacts and additional multi-particle states that could be created by the lattice 
fermion field operator \1/ beyond the single quark state. Writing 

-G( Po ,p=0) = , 1 (A2) 
a t %<y A + B 

and using the complex variable 

z = e iatpo = e~ atE , (A3) 

one finds 



17 



One can show that the pole, 



dz 



-ij A + B 



2mz (A + iB)(A — iB) 



(A4) 



Z\ = e 



-a* Mi 



(A5) 



corresponding to a physical positive energy particle obeys (A — iB)\ z=Zl = or B\ z=Zl = 
—iA\ z=Zl . The contribution to G(t, 0) from the residue at this physical pole is then given 
by, 



?* /at H7o^ + £) 



(70 + 1) -Aft* 

2 



— z 



2 (A' - iB') 



(A6) 



Using 



dz 



df 



d(a t po) 



p =iMi 



and comparing with eq . (|AT|) , one finds, 



(A7) 



d(atpo) 



(A - iB) 



Po =iM 1 



(A8) 



For <Sq234°^ one has at tree-level, 



A (aniso) = sin ( am ) 

= a t m + x-Xcos(a tPo ) 



(A9) 
(A10) 
(All) 



and 



-a t Mi 



(a t m + x) ~ \j (a t m + x) 2 + 1 - X 2 



or equivalently, 



(A12) 



,a t Mi 



(a t m + x) + V ( a * m + X) 2 + 1 - X 2 
X + l ' 



(A13) 



As explained in reference |L8| m = tjiq — m c and m c is the value of mo that gives a massless 
pion. From eq. (|A8|) one finds, 



y(0),aniso 

q 



1 



coshiatM^ + xsinhiatMi) ^(a t m) 2 + 2(a t m)x + I 



(A14) 
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For the isotropic action with higher time derivatives the formulas are more compli- 

cated. 

A {iso) = ^sin(a t po) - \sin{2a t p ) (A15) 
3 6 

B (iso) = a t m + -(1 - cos(atpo)) - \sin 2 {a tPo ). (A16) 
3 6 



(A17) 



Eq(P) leads to 



|e atMl - |(cos/i 2 (aiM!) + sinh 2 (a t M l ) + sinh(a t Mi)cosh(a t M 1 )) 

The tree-level physical pole = e~ atMl is the solution to, 

z 4 - (24a t m + 30)z 2 + 32z - 3 = (A18) 

that evolves smoothly from Z\ = 1 at a t m = 0. A lengthy closed expression for £i as 
a function of c^m can be obtained (using for instance Mathematica) , however we do not 
consider it worthwhile reproducing it here. It is easier to plug in specific values for a t m into 
eq. ([A18|) before solving for the roots. At small a t m the physical pole is well approximated 

by 

l = e «t^ = l + ( atm ) + I( atm )3_ A(a t m) 4 + 0((a t mf). (A19) 
Z\ 6 24 

For the mass parameter values used in the current simulations, the sum of just the first three 
terms in ( A19| ) differs from the exact solution to eq.( |Al£| ) by only 1%. 

In order to get explicit values for Zfp, one needs to know a t m = a t (m — m c ). Since we 
have data at only one light quark mass, we do not have a nonperturbative estimate for atm c 
based on a vanishing pion mass. We approximate atm c using perturbation theory and find 
ITS]] , a t m fa 0.635 and a t m ~ 0.196 respectively for the isotropic and anisotropic lattices. 
This leads to 



7 (0),aniso 

which are the values used in section IV. 



z (0),iso = ^^_ (A2Q) 

9 1.2226 
1 



9 1.3286 v ; 
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TABLES 



TABLE I. Simulation Details. 



isotropic 



anisotropic 



lattice size 
# configs 

P 

Landau link uq 
Xo 

X = a s /a t 
Co 



a t m 
P/V 

typical # of BiCGstab iters 
a s M 



3 3 x 20 
200 
1.719 



0.797 p| 



1 
1 

0.8(1) GeV 
0.8(1) GeV 
1.15 
0.725(5) 
140 - 160 
6.5 



S 3 x 48 
200 
1.8 



u s =0.721 « t =0.992 |f| 
6.0 



5.3 |7] 
0.82 
0.7(1) GeV 
3.7(4) GeV 
0.39 
0.726(6) 
220 - 270 
7.0 and 2.0 



TABLE II. Some Spectrum Results. The first errors are statistical and the second are estimates 
for errors due to uncertainties in the scale. The hyperfine splittings have not been adjusted for 
incomplete tuning of the heavy quark mass (see text for adjusted numbers). 





isotropic 


anisotropic 


Experiment 


light-light 




"pion" mass 


0.856(3)(86) GeV 


0.840(7)(84) GeV 




"rho" mass 


1.179(10)(118) GeV 


1.158(7)(116) GeV 




heavy-light 




Bs 


5.65(31)(57) GeV 


4.80(20)(48) GeV 


5.369 GeV 


D s 




2.04(5) (20) GeV 


1.969 GeV 


B* s -B s 


46(3) (5) MeV 


45(2) (5) MeV 


47.0(26) MeV 


D* s ~ D s 




143(4) (14) MeV 


143.8(4) MeV 
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FIG. 1. Fractional errors in B s meson correlators for momentum (0,0,0) and (1,1,1). 
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FIG. 2. Fractional errors in B s meson correlators for momentum (0,0,2) and (0,0,3). 
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FIG. 4. Effective masses for the (1,1,1) momentum B s meson correlators in (i) lattice and (ii) 
physical units. The isotropic points in (ii) have been shifted down for clarity. 
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FIG. 7. Effective masses for the (1,1,1) momentum pion correlators in (i) lattice and (ii) 
physical units. The isotropic points in (ii) have been shifted down for clarity. 
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FIG. 9. Effective masses for the (1,1,1) momentum rho correlators in (i) lattice and (ii) 
physical units. The isotropic points in (ii) have been shifted down for clarity. 
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FIG. 25. T^>{t) of eq.pl) for pion momentum (1,1,1) versus time in lattice units. 
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FIG. 26. T^\t) of eq.© for pion momentum (0,0,2) versus time in lattice units. 
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FIG. 28. The form factors fo(q 2 ) and f+(q 2 ) for B s decays from isotropic lattices. See Table 
II for actual values of the decaying heavy meson and daughter meson masses. 
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FIG. 29. The form factors fo{q 2 ) and f+(q 2 ) for B s decays from anisotropic lattices. See Table 
II for actual values of the decaying heavy meson and daughter meson masses. Some points have 
been shifted horizontally for clarity. 
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FIG. 30. The form factors fo(q 2 ) and f+(q 2 ) for D s decays from anisotropic lattices. The 
"burst" shows an experimentally determined value of f+(q 2 = 0) for the decay D° — ► K~l + v [34]. 
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